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Non-Abelian statistics is a phenomenon of topologically protected non-AbeUan Berry phases as 
we exchange quasiparticle excitations. In this paper, we construct a Zn rotor model that realizes 
a self-dual Zn Abelian gauge theory. We find that lattice dislocation defects in the model produce 
topologically protected degeneracy. Even though dislocations are not quasiparticle excitations, they 
resemble non-Abelian anyons with quantum dimension a/]V- Exchanging dislocations can produces 
topologically protected projective non-Abelian Berry phases. The dislocations, as projective non- 
Abelian anyons can be viewed as a generalization of the Majorana zero modes. 
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Introduction — Searching for Majorana fermions (or 
more precisely, Majorana zero modes) in condensed mat- 
ter systems have attracted increasing research inter- 
ests recently. [I]-[9] But what really is the Majorana zero 
mode? In fact, the so called "Majorana zero mode" is 
actually a phenomenon of topologically protected degen- 
eracy in the presence of certain topological defects (such 
as vortices in 2D + ipy superconductors [2| i3j). In 
the race for finding Majorana zero modes, much atten- 
tion has been paid to the fermion systems. [4]-[9] However 
the boson/spin systems also have topologically protected 
degeneracies, p^HTO] which may also be ascribed to Ma- 
jorana zero modes or their generalizations. 

An ID example of emergent Majorana zero modes in 
the spin system arises from the transverse field Ising 
chain, [20l [21] whose ground state degeneracy in the fer- 
romagnetic phase can be viewed as the Majorana zero 
modes at both ends of the chain. A 2D example is found 
in the toric code model [151 Hi]? where lattice disloca- 
tions are braided and fused as if they were Majorana zero 
modes [22l[23| which resemble non-Abelian anvons ^24] - [26] 
of quantum dimension a/2. The toric code model can be 
generalized to a Z^v rotor model, whose low energy ef- 
fective theory is a self-dual Zn gauge theory. [27'/2 8 In 
this paper, we study the topologically protected degen- 
eracy associated with the extrinsic topological defects, 
namely lattice dislocations in the Zn rotor model, and 
found that these defects are of quantum dimension a/TV, 
which can be viewed as a generalization of the "Majorana 
zero mode" . Braiding topological defects with protected 
degeneracy will lead to topologically protected projective 
non-Abelian Berry phase, which may allow us to perform 
decoherence free quantum computations. ^26^ 

We like to remark that the dislocations in our Zn rotor 
model are not non-Abelian anyons, since the non-Abelian 
anyons must be excitations of the Hamiltonian, while the 
dislocations are not the excitations in this sense. The dis- 
locations do not really carry non-Abelian statistics since 
the non-Abelian Berry phase from exchanging disloca- 



tions is topologically protected only up to a total phase. 
We say the dislocations carry a projective non-Abelian 
statistics. [29, 30 An other example of projective non- 
Abelian statistics for dislocations in fractional quantum 
Hall states on lattice can be found in Ref. 31, 

Ztv plaquette model — The Z^v plaquette model is a ro- 
tor model on a two-dimensional square lattice (see Fig.[T]). 
On each site define a Z^ rotor with N basis states |m^), 
labeled by the angular momentum rui = 0, 1, • • • , (N — l). 
For each rotor, introduce Ui to measure the angular mo- 



mentum Ui\mi) 



\mi) with On = 2ii/N^ and 



Vi to lower the angular momentum by one Vi\mi) = 
\{mi — l)mod n)- Both Ui and Vi are unitary operators 
UjUi = V^Vi = 1, satisfying ViUi^ = e'^^^^^' Ui'Vi. 

The Zat plaquette model is given by the Hamiltonian 



H = - ^Op^ h.c, 



(1) 



where the operator Op describes a kind of ring coupling 
among the rotors on the corner sites of each plaquette p. 



Op = ITp) = 



(2) 



Here we adopt the graphical representation for the oper- 
ators: Ui =\, Vi = /, U\ = \, V} = by drawing di- 
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FIG. 1: (Color on line.) Even x even periodic lattice with 
plaquettes colored in a check board pattern: red and darker 
(blue and lighter) plaquette will be called even (odd). Each 
directed string represent a product of Ui and/or Vi operators 
on the sites along the string. The operator on each site is 
specified by the string direction (see text). 
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rected strings going through the site. Because these oper- 
ators only connect diagonal plaquettes, a string starting 
from the even plaquette will never enter the odd plaque- 
tte (and vice versa). So we can locally distinguish two 
different types of strings: e-string (m-string) if it lives in 
the even (odd) plaquettes (see Fig.jl]). The assignment 
of even/odd to the plaquettes can be reversed under the 
translation of one lattice spacing, so the interchange of e- 
and m-strings could be realized by the curvature of the 
lattice as will be seen later. 

The Zat plaquette model Eq. ([T]) is exact solvable, as 
evidenced from the commutation relation [0^,0^/] = 0, 
as OpOp> = @J = e^^^e-*^^g@ = 0^.0^ for ad- 
jacent p and p' ^ where the overlay of strings indicates 
the ordering of the operators, such as ViUi = X 
UiVi = X, with the algebra >2 = e^^^)^. 

Every Op operator has N distinct eigenvalues e^^^^p 
labeled by qp = 0, 1, • • • , (A^ — 1), as inferred from the 
fact that = 1. qp denotes the (generalized) Ztv 
charge hosted by the plaquette p. If the plaquette is 
even (odd), we may call it e-charge (m-charge). The en- 
ergy will be minimized if all Op's take the eigenvalue 1 
{% = 0)- Therefore the ground states are the common 
eigenstates that satisfy Op|grnd) = |grnd) for ally's, and 
is free of any charges. 

Intrinsic anyon excitations — The excited states can 
be obtained by applying open string operators to the 
ground state, which create opposite Zn charge excita- 
tions in pairs at both ends of the string. These exci- 
tations and can be detected by the close string opera- 
tor (like Op) surrounding them in the counterclockwise 
direction. Take S in Figjl] for example, Op5'|grnd) = 
1^^- • • Igrnd) e*^^ • • -l^jgrnd) = e^^^^lgrnd), show- 
ing that a charge qp = +1 is created at the end of the 
open string by the action of S. One can show that the 
opposite charge = — 1 is created at the other end. 

Because Zn charge excitations are the ends of open 
strings, their statistics are inherited from the algebra of 
the string operators. According to = e*^^^, braiding 
a qe e-charge with a qm m-charge would acquire a phase 
^iieN/'2)qeqm ^ In this sense, these excitations are Abelian 
anyons. However we must stress that these anyons are in- 
trinsic^ as they are collective motions of rotors, described 
by the excited state within the rotor Hilbert space. This 
is to be distinguished from the extrinsic anyons intro- 
duced later as lattice dislocations, which does not be- 
longs to the rotor Hilbert space. Note that both the 
phase e*^^^/^^*^^*^"^ and the excitation energy are invari- 
ant under the exchange of e and m. This manifests the 
self-duality of the Zn plaquette model, and can be real- 
ized by lattice translation. 

Ground state degeneracy — The degeneracy of the 
ground states of the Z^ plaquette model depends on 
the topology of the lattice. Let us consider the torus 
topology by setting the model on a L^^ x Ly sized lat- 



tice with periodic boundary condition in both directions. 
The total number of states is A^^^^te^ with A^site = L^Ly 
being the number of sites. To count the ground states, 
we note that they are constrained by Vp : Op = 1. Con- 
sider a particular Op operator and the subspaces labeled 
by its different eigenvalues. Those subspaces all have 
the same dimension, because any open string operator 
that ends in the plaquette p can be used to perform 
a unitary transform that rotates these subspaces into 
each other. So each time imposing Op = 1 on a par- 
ticular plaquette will reduce the available Hilbert space 
dimension by a factor of N. However the Op operators 
are not independent. Because e-charges (m-charges) are 
created in opposite pairs, summing over the lattice, e- 
charges and m-charges must be neutralized respectively, 
i-e. Upeeven^P = IIpGodd = 1- This is truc on an 
even X even lattice (i.e. both Lx, Ly are even), which 
reduces the number of independent Op constrains to 
(A^piaq — 2), with TVpiaq = LxLy being the number of 
plaquettes. So after restricting the full Hilbert space to 
the ground state subspace, the remaining dimension is 
_ ^^2^ meaning the ground state degener- 
acy of the Zn plaquette model is N'^ on the even x even 
lattice. However for the even x odd or odd x odd lattices 
(i.e. Lx or Ly is odd), e-string and m-string can be con- 
tinued into each other by going along the odd direction, 
thus e-charge and m-charge are made identical. So they 
are no longer required to be neutralized respectively, but 
only neutralized as a whole. Therefore we only have one 
relation flp = 1^ which reduces the number of inde- 
pendent Op constraints to (A/'piaq — 1), and the resulting 
ground state degeneracy will be 7V^^i*^~^pi^^+^ = N. 

To summarize, the ground state degeneracy of the Zn 
plaquette model on a torus follows from the general for- 
mula 

GSD = A/'7V^^^^^-^pl-^ (3) 

where JV denotes the number of species of the intrin- 
sic excitations that are supported by the lattice topol- 
ogy. On the even X even lattice, we have totally JV = N'^ 
distinct excitations by combination of e- and m-charges. 
When it comes to the even x odd or odd x odd lattice, e- 
and m-charges are no longer distinct, and the number of 
excitation species is reduced io M = N . The topologi- 
cal order in the ground state is now evidenced from the 
protected ground state degeneracy on torus, ^Olln] and 
from the dependence of the ground state degeneracy on 
the parity of the lattice periodicity. 

Dislocations — One can change the lattice periodicity 
by first generating a pair of edge dislocations with op- 
posite unit length Burger's vectors, and moving them in 
the direction perpendicular to their Burger's vectors all 
the way around the lattice, then annihilating them as 
they meet again at the periodic boundary. During this 
process, the ground state degeneracy must have changed. 
This motivates us to introduce dislocations as shown in 





FIG. 2: (Color on line.) (a) Lattice with a pair of disloca- 
tions, marked out by H and h. Plaquettes on the branch cut 
are colored by violet. Periodic boundary conditions are as- 
sumed in both direction by sticking the dashed edges with 
the solid edges on the opposite side. The ring operators Op 
are redefined around the pentagonal plaquette. C operators 



denote large close strings looping around the lattice, (b) Pla- 
quette to site mapping. The site that is not mapped to is 
marked by a black dot. 



FIG. 3: (Color on line.) Diffeomorphism of the torus with 
a pair of dislocations at c and d. Expand the branch cut 
(violet dashed line) between c and d into a hole, with two 
edges marked by a and /3. Separate the e- and m- layers. 
Unwrap both layers by cutting along the large loops around 
the torus. Rotate one layer to glue the /3 edges together along 
the marked direction. Glue the other edges and rewrap into 
a double torus. 



Fig. |2(a)| to probe the topological order by looking at the 
degeneracy associated to them. With dislocations, one 
can no longer globally color the plaquettes consistently. 
Branch cuts must be left behind between each pairs of 
dislocations. Going around a dislocation exchange the 
e- and m-charges, as e- and m-strings are transmuted 
into each other across the branch cut. The self-duality is 
made explicit by dislocations. 

In the presence of dislocations, the Z^r plaquette model 
is still defined by the Hamiltonian in Eq. ([T]), with the 
same ring operator Op in Eq. ([2| for quadrangular plaque- 
ttes (including those on the branch cuts). Only around 
the pentagonal plaquettes (at the dislocations), the ring 
operator Op should be redefined as 



Op = 



4S = -e'^UiV2UlV^U,Vl (4) 



The phase factor —e^^^l^ is to guarantee that O^ = 1 
holds for the pentagonal plaquette as well. The pentag- 
onal ring operator Op commutes with all the other ring 
operators, so the exact solvability of the model is pre- 
served. The ground states are again common eigenstates 
of \/p : Opigrnd) = |grnd). The dislocations are topo- 
logical defects that do not belong to the model Hilbert 
space. To distinguish from those intrinsic Zn charges, 
we will call the dislocations as the extrinsic defects. 

With the branch cuts, e-charge and m-charge are indis- 
tinguishable, so the species of intrinsic excitations count 
to Af = N. According to Eq. (|3|, the ground state de- 
generacy will be given by 7V^site-^piaq+i [^i general. To 
count the number of sites and plaquettes, we first es- 
tablish a correspondence between them by mapping each 
plaquette to its bottom-left corner site, as indicated by 
the arrows in Fig. 2(b) Between a pair of dislocations. 



correspondence (see Fig. |2(b)] ), so the introduction of ev- 
ery pair of dislocations will give rise to one extra site 
(with respect to the number of plaquettes). Therefore 
if there are n dislocations on the lattice, there will be 
A/'site — A'piaq = n/2 morc sites than plaquettes, and the 
ground state degeneracy of the Zat plaquette model will 
be GSD = 7V^/2+i. 

This ground state degeneracy is topologically pro- 
tected indeed. To better understand the topology, we 
start from the even x even periodic lattice without dislo- 
cations, i.e. a torus with no branch cut. In this case, 
the e-strings and m-strings are distinct, and can never 
be deformed into each other, as if they were living on 
two different layers of the torus. So the topological space 
is the disjoint union of two separate torus. Introducing 
a pair of dislocations, the two layers will be connected: 
strings on one layer can be carried on into the other layer 
through the branch cut. So the topological space be- 
comes a doubled torus under the diffeomorphism [32j as 
shown in Fig. [3) 

All the operators that act within the ground state sub- 
space are closed-string (cycle) operators, as they com- 
mute with the Hamiltonian. Note that the contractable 
cycles act trivially (as Op = 1). Only non-contract able 
cycles can be used to label the different ground states 
and to perform unitary transforms among them. On 
the double torus topology as in Fig. |4(a)[ one can spec- 
ify 4 non-contractable cycles: Cex^ Cey^ Cmx^ Cmy^ as 
the canonical homology basis. Their operator forms are 
given explicitly according to their graphical representa- 
tions depicted in Fig. 2(a) We now study the represen- 



tation of these cycle operators in the ground state sub- 
space. First we find the following commutation relations 



and two independent algebras CeyC^ 



ey^mx 



^mx^eyi 

CexCmy Each algebra requires an A"- 
only one of them will hold a site that has no plaquette dimensional representation space, so the 4 cycle operators 



CmyCe, 
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FIG. 4: (C olor on line.) Diffeomorphism of strin g operators 
on the torus (a) with a pair of dislocations, or (b) with 3 pairs 
of dislocations. 



together requires A^^ -dimensional representation space, 
which must have completed the ground state subspace, 
since all the non-contractable cycles can be generated by 
these 4 basis cycles. Therefore the ground states are A^^- 
fold degenerated, and each of them corresponds to a ba- 
sis in the representation space. Any perturbation of the 
Hamiltonian that is non-zero only in a compact region 
will not change the ground state degeneracy, since the 
non-contractable cycle operators that avoid the compact 
region still commute with the Hamiltonian. 

The above can be generalized to the case with any 
number of dislocations. Consider n dislocations with 
n/2 branch cuts. Following the similar cut-and-glue pro- 
cedures in Fig. [3) the topological space will be a genus 
g = n/2 + 1 surface as in Fig. |4(b)| on which one can 
choose g pairs of non-contractable cycle operators Ca and 
(a = I,-- - ,^), such that [C,, C5] = K,^] = 
and CaC^ = e^^^^^^C'^Ca- These operators spans a 
A^^-dimensional representation space isomorphic to the 
ground state subspace. Therefore the ground state de- 
generacy of the Zat plaquette model with n dislocations 
is GSD = A^^/2+^, which is consistent with our previous 
result. Each dislocation contributes to the ground state 
degeneracy by a factor of \fN . Thus the dislocations re- 
semble non-Abelian anyons of quantum dimension a/TV, 
as described in Ref. 32 . Braiding the dislocations leads 
to topologically protected projective non-Abelian Berry 
phases. We see that projective non-Abelian anyon can 
emerge from an Abelian model as the extrinsic topologi- 
cal defects, such as lattice dislocations. Those projective 
non-Abelian anyon can be used to perform topological 
quantum computations, [26 but not universally since the 
square of the quantum dimension is an integer. [33 

Parton approach — For the N = 2 case, the quantum 
dimension \/2 implies that the extrinsic anyons are Ma- 
jor ana fermions. To expose the Major ana fermion ex- 



plicitly, we evoke the parton projective construction, in 
which 4 Majorana fermions rff (a = 1,2,3,4) are in- 
troduced on each site i, obeying the ant i- commutation 
relation {rjf^r]^} = (^ij^Q;/3.^3 Under the constraint 
1/4, the rotor operators can be expressed 



as Ui = iT]lr]f^ Vi = irjfrjf. Then the Z2 plaquette 
model can be mapped to an interacting fermion model, 
which has a "mean-field" description given by i^mean = 
— + ^'^') with the ansatz Sij = ±1 on each 

bound, where A^^^+^ = ir]jr]f_^^ and A- • - — '^^'^^^ 



Let be a free fermion ground state of i^mean, and 

V = Y\i\{^ ^ "^ViViViVi) be the projection operator to 
the physical Hilbert space of rotors. All the eigenstates 
of the Zat plaquette model can be obtained by the pro- 
jective construction as V\{sij}). To obtain the ground 
states, {sij} must satisfy the fiux configuration given by 
Op = 1, which has totally 4 gauge inequivalent solutions 
on a torus. Given a particular {s^j}, all the Majorana 
fermions will be paired up across the bound, except for 
the dangling Majorana fermion at the dislocation site. 
If there are n dislocations in the system, there will be 
n dangling Majorana zero modes, which leads to a 2^/^ 
fold degeneracy in the free fermion ground states. So al- 
together we have 4 x 2"^/^ fermion states to be projected 
from, half of which will be projected to nothing due to 
their odd fermion parity. Therefore the resulting physical 
ground states add up to 4 x 2"^/^/2 = 2^/^+^, consistent 
with our previous formula. The above discussion has 
shown that the a/2 quantum dimension of the extrinsic 
anyon actually originated from the dangling Majorana 
fermion, or the Majorana zero mode, at the dislocation 
site. It has been shown that exchanging Majorana zero 
modes will lead to non-Abelian Berry phase, which sup- 
ports our conjecture that exchanging dislocations in our 
Ztv plaquette model leads to protected (projective) non- 
Abelian Berry phase. 

In conclusion, we studied the phenomenon of topolog- 
ically protected degeneracy and topologically protected 
projective non-Abelian Berry phases produced by extrin- 
sic topological defects (such as dislocations) in a Z^v ro- 
tor model. We find that these dislocations are projective 
non-Abelian anyons with quantum dimension ^/N. For 
A/" = 2, such a result can be re-derived from a parton 
construction where the dislocations can be identified as 
Majorana zero modes. For higher N {N > 2), the pro- 
jective non-Abelian anyons {i.e. the dislocations) can be 
viewed as a generalization of the Majorana zero modes. 
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